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Abstract 



A theoretical prediction is given for the spin and orbital angular momentum distribution 
functions of the nucleon within the framework of an effective quark model of QCD, i.e. the 
chiral quark soliton model. An outstanding feature of the model is that it predicts fairly small 
quark spin fraction of the nucleon AS ~ 0.35, which in turn dictates that the remaining 65 % 
of the nucleon spin is carried by the orbital angular momentum of quarks and antiquarks at 
the model energy scale of Q 2 ~ 0.3 GeV 2 . This large orbital angular momentum necessarily 
affects the scenario of scale dependence of the nucleon spin contents in a drastic way. 
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1 Introduction 



The recent renewal of interest in polarized lepton-hadron and hadron-hadron deep-inelastic 
scatterings greatly owes to the so-called "nucleon spin crisis" aroused by the EMC experiment 
|U]. According to the latest SMC analysis @, the intrinsic quark spin carries only about 20 % of 
the total nucleon spin. The widely-believed resolution of this puzzle is based on an assumption 
of large gluon polarization ||. In fact, the recent next-to-leading-order (NLO) analyses of the 
polarized deep-inelastic-scattering (DIS) data for gi(x,Q 2 ) appear to indicate a positive value 
for the gluon polarization [|J . Nevertheless, the situation is far from conclusive, since the result 
of these analyses is rather sensitive to the input shapes of the polarized distributions as well 
as the uncertainties of the available polarized DIS data. 

Although somewhat nonstandard, there is another scenario that also leads to qualitative 
resolution of the nucleon spin problem. This scenario was first introduced on the basis of the 
Skyrme model ||. Later, the idea was developed into a more quantitative one within the 
framework of the chiral quark soliton model (CQSM) pi, [71]. The CQSM is an effective quark 
model of QCD, which was first introduced by Diakonov and Petrov, based on the instanton- 
liquid picture of the QCD vacuum f|,[f|. Being in conformity with Witten's idea of large 
N c QCD JTI]| , it shares many features in common with the Skyrme model ||12||,||13|| at least in 
the ideal limit of N c = oo. (For review, see flOf .) At the subleading order of 1/N C expansion, 
however, several crucial differences have been found to exist between them p)| , ||14|| , |I5| . Among 
others, most important would be the following two. The first is the 1/N C correction (or the first 
order rotational correction in the collective angular velocity Q) to the isovector axial-vector 
coupling constant of the nucleon, which definitely exists in the CQSM but is entirely missing 
in the Skyrme model Jl5[ . (In our opinion, this solves the so-called 11 g a problem" in the latter 
nj.) The second concerns the subject of our central interest here, i.e. the predictions for 



the quark spin fraction AS of the nucleon. As was first noticed by Brodsky et al. ||, AE 
identically vanishes in the naive Skyrme model. At first sight, this appeared to be a favorable 
feature of the model since it looks consistent with the result of the EMC measurement lUl . As 



we have argued repeatedly however, this outstanding prediction of the Skyrme model 



should be taken with care, because it just arises from the fact that this model cannot correctly 
describe the physics at the subleading order of 1/N C expansion. In fact, the presence of the 
next-to-leading order contribution to AE is again a distinguishable feature of the CQSM as 
compared with the naive Skyrme model Nonetheless, qualitative similarity between the 
two models still survives and the prediction of the CQSM for AE remains to be fairly small, 
i.e. AE ~ 0.35, as compared with other effective models of the nucleon like the nonrelativistic 



quark model or the MIT bag model. What carries the rest of the nucleon spin ? Within the 
CQSM, the answer to this question is quite simple. It must be the orbital angular momentum 
of quarks and/or antiquarks, since it is an effective theory with the effective quark degrees of 
freedom only Q]. 

A natural question is whether the effective quark degrees of freedom are enough to describe 
the relevant physics. One might suspect that the gluon polarization or its orbital angular mo- 
mentum may play important roles even at the low energy hadronic scale. There have been 
some attempts to introduce the explicit gluon degrees of freedom into the CQSM in reference 
to the underlying instanton-liquid picture of the QCD vacuum ||17||,||18|. Unfortunately, the 
situation is far from conclusive yet. Here we simply assumes that the explicit gluonic degrees 
of freedom are not crucial at least for the polarized distribution functions at the low renormal- 
ization point, which is taken to be the starting energy of the DGLAP evolution equation. That 
this is not a meaningless assumption can, for instance, be convinced from our recent analysis 
of the longitudinally polarized distribution functions of the nucleon and the deuteron in com- 
parison with the recent EMC and SMC data , |20|] . This analysis starts with the predictions 
of the CQSM for the longitudinally polarized distribution functions for quarks and antiquarks, 
which are taken to be input distributions given at the model energy scale of Qf ni = 0.3 GeV. 
The polarized distribution functions for the gluon are simply assumed to be zero at this initial 
energy scale. After solving the standard evolution equation at the NLO, the predictions of 
the model were then compared with the recent EMC and SMC data for gi(x, Q 2 ), g™ (x, Q 2 ), 
and gf(x,Q 2 ) given at Q 2 = 5 GeV 2 0. Despite the rather drastic assumption of zero gluon 
polarization at the hadronic scale, the theory was shown to reproduce all the qualitatively 
noticeable features of the experimental data. 

Encouraged by this success, now we push forward our attempts to clarify the spin contents of 
the nucleon. That is, the main purpose of the present study is to give a theoretical prediction for 
the orbital angular momentum distribution functions for quarks and antiquarks in the nucleon. 
After combining those with the previously-obtained quark spin distribution functions, we can 
investigate the scale dependence of the full nucleon spin contents by using the recently-derived 
evolution equation at the leading order (LO) [p^j , ||23|| . (We however recall fact that there is a 
criticism against this LO evolution equation for the orbital angular momentum distributions 
J33J.) Although our analysis is bound to one basic assumption that the gluon polarization as 
well as the gluon orbital angular momentum play no significant role at least at the hadronic 
scale of Q 2 < 0.3 GeV, it is expected to provide us with quite a unique scenario for the scale 
dependence of the nucleon spin content, while keeping satisfactory agreement with the available 
high-energy data for g\ (x, Q 2 ),g™(x, Q 2 ), and gf (x, Q 2 ). 
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We shall explain in the next section the theoretical framework which our calculation of 
the orbital angular momentum distribution functions for quarks and antiquarks are based on. 
Next, we discuss the scale dependence of the nucleon spin contents including the orbital angular 
momentum based on the given input distributions at the low renormalization point. Finally, 
we summarize our findings in the last section. 



2 The theoretical framework 

How to treat the orbital angular momentum has been one of the most controversial theoret- 
ical problems in hadron spin physics. It is well known that the most natural definition of 



angular momentum is not gauge invariant |25 ]. Recently, Bashinsky and Jaffe proposed a 



gauge-invariant formulation of the angular momentum for quarks and gluons which reduces 
to natural definition in a particular gauge, i.e. the light-cone gauge A + (x) = [26]. It was 



also suggested |27f that this gauge invariant formulation may give an extended support to 



the recently derived LO evolution equation for the orbital angular momentum distributions in 
the light-cone gauge f2~2"f,||23||. According to the authors of [24], however, nonlocal operators 



with dependence on spatial coordinates have not been seen in factorization of hard forward 
scattering processes, and in particular inclusive deep-inelastic scatterings do not depend on 
such types of operators. Accodingly, the evolution of the orbital angular momentum distribu- 
tions would become extremely complicated in the light-cone gauge contrary to the results of 



22 1, [23]. On the other hand, there is another definition of orbital angular momentum distribu- 



tions, which maintains gauge invariance at the cost of losing full Lorentz covariance P3 |. This 
latter definition has its meaning only in a limited class of coordinates in which the nucleon has 
a definite helicity. An advantage of the second definition is that it can be extracted from the 
polarized and unpolarized quark distributions and the off- forward distributions E{x) so that it 
can in principle be measured. What is not fully clarified yet is the relation between these two 
definitions of angular momentum distributions. In sum, the proper definition of the orbital 
angular momentum distributions as well as the corresponding evolution equation in QCD is 
still an unsettled issue being subject to a debate. One should therefore take the following 
investigation of the orbital angular momentum distribution within an effective theory of QCD 
bearing all these affairs in mind. 

Since the CQSM at the present level of approximation contains no gluonic degrees of 
freedom at least explicitly, it seems natural to start with the following naive definition of the 
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quark orbital angular momentum distribution function : 

I POO rj\ 

q L (x) = -=— / — e lXx <P^|^(0)^(^ 2 -x 2 9 1 )^ + (An)|P5 3 >, (1) 



y/2p+ J-oo 2 IT 

Here p^ and n M are two light-like (null) vectors, having the properties 

p~ = 0, n + = 0, p 2 = n 2 = 0, p ■ n — 1, (2) 
while ■0+ is a component of the quark field ip defined through the decomposition 

= (P+ + P-)V = 0+ + 0-, (3) 

by the projection operators P± = |7 =F 7 ± with 7 ± = (l/y/2) (7 ±7 3 ). Extending the definition 
of distribution function to interval —1 < x < 1, the relevant antiquark distributions are 

given as 

= (0 < a; < 1). (4) 

Naturally, from more general viewpoint of underlying color gauge theory, i.e. QCD, the naive 
definition (1) is not gauge invariant and holds only in the light-cone gauge. As was already 
explained, this definition of the quark orbital angular momentum is not free from problems, 
especially when we are to use it as a initial distribution of scale evolution. We expect that 
this is not so serious for our investigation below, since our main concern there is the qualita- 
tive behavior of the quark orbital angular momentum distribution at the model scale, which 
naturally has a strong correlation with the quark spin distribution at the same scale. On the 
other hand, the study of its scale dependence should be taken as semi-qualitative one. 

The basis of our theoretical analysis is the following path integral representation of the 



nucleon matrix element of bilocal quark operator p8 



< N(P) I + (O) Oip(z) I N(P) > = — J ' d 3 x d 3 y e~ iP x e iP y J Vir J Vip 

x j N (^, x ).^\0)O^(z)-4(~,y)e i f d3 * c °^, (5) 

where 

Ccqm = j>(i?-Mj-* T -W*)1>, (6) 
is the basic lagrangian of the CQSM , while 

Ms) = e 0tt, " ,ow - r &,TT 3 c VWiW • "^.W. (7) 
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is a composite operator carrying the quantum numbers JJ 3 , TT 3 (spin, isospin) of the nucleon, 

f ... f 

where a, is the color index, while ^ c is a symmetric matrix in spin-flavor indices /j. We 
start with a stationary pion field configuration of hedgehog shape : 

n(x) = UrF(r). (8) 

Next we carry out a path integral over tt(x) in a saddle point approximation by taking care 
of two zero-energy modes, i.e. the translational zero-modes and rotational zero-modes. Under 
the assumption of "slow rotation" as compared with the intrinsic quark motion, the answer 
can be obtained in a perturbative series in Q, which can also be regarded as a 1/N C expansion. 
The first nonvanishing contribution to qi{x) arises at the O^ 1 ) term of this expansion, since 
the leading 0(Q°) term vanishes identically due to the hedgehog symmetry. According to the 
general formalism derived in the answer is given in the following form : 



Aq L (x) = Aq[ A > B \x) + AqZ(x), (9) 

where 

Aq { L AB} = <J 3>pr .M N ^- F \ F <n\T 3 \m><m\{l + 1 °^)L 3 5 n \n> 

1 m=all,n<0 ^ m ^ n 

= - < J3 > P T ' M n ~r Y) T-i 1 T-i <n\r 3 \m><m\(l + 7V )L 3 8 n \ n >, (10) 

i a — o £j m — &n 



m=alLn>0 



and 



Aq c L (x) = < J 3 > pT -j- J2 <n\r 3 (l + 1 ° 1 3 )L 3 5 n \n> 

aX A1 n<0 

= - <^3> PT Y,< n \T3(l + l°l 3 )L 3 5 n \n>, (11) 

aX A1 n>0 

with L3 = — i (a; 1 ^ 2 — x 2 d l ) and 5 n = 5{xMn — E n — p 3 ), while (O) ^ denotes a matrix element 
of a collective space operator O with respect to the proton in the spin up state along the z-axis, 
i.e. 

<0>rf = / ^){y^ ^^ = (p,S 3 = l/2\0\p,S 3 =l/2). (12) 

Here, | m > and E m are the eigenstates and the associated eigenenergies of the static hamil- 
tonian H = — ia • V + (3M e* l5 ' T ' r F ^ with the hedgehog background. Note that both 
of Aqjf' B ^ (x) and q < l{x) are represented in two equivalent forms, i.e. in the occupied form 
(n < 0) convenient for the numerical calculation for x > and in the nonoccupied form (n > 0) 
convenient for x < 0. 
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In the actual numerical calculation, the expression of any physical quantity is divided into 
two pieces, i.e. the contribution of what we call the valence quark level ( it is the lowest energy 
eigenstate of the static hamiltonian H, which emerges from the positive energy continuum) 
and that of the Dirac sea quarks (the latter is also called the vacuum polarization contribution) 
as explained in Since the latter contains ultraviolet divergences, it must be regularized. 
Here we use the so-called Pauli-Villars regularization scheme p8 |. The regulator mass Mpy 
of this cutoff scheme is determined uniquely from the physical requirement that the effective 
meson action derived from (6) gives correct normalization for the pion kinetic term. Using 
the value of M — 375 MeV, which is favored from the phenomenology of nucleon low energy 
observables, this condition gives Mpy ~ 562 MeV. Since we are to use these values of M 
and Mpy, there is no free parameter additionally introduced in the calculation of distribution 
functions. 



3 Numerical Results and Discussion 

We first show in Fig.l the theoretical predictions for the orbital angular momentum distribution 
functions of quarks and antiquarks as well as the isosinglet quark polarization. (The latter was 
already given in [jTSfl , but it is shown here again for emphasizing how the x-dependences of the 
two distributions are different.) Here the distribution with negative x are to be interpreted as 
the antiquark distribution according to the rule (4) and the similar relation for Au(x) + Ad(x). 
The long-dashed and dash-dotted curves respectively stand for the contribution of the discrete 
valence level and that of the negative energy Dirac sea, while their sum is shown by the 
solid curves. As already noticed in |TI| , the contribution of Dirac sea is not so significant for 



Au(x) -\-Ad(x). The situation is quite different for the orbital angular momentum distribution 
qh{x). One clearly sees that the contribution of the polarized Dirac vacuum being peaked 
around x ~ is dominating over that of the discrete valence level. This can also be convinced 
from the first moment defined by 

L q = J ^q L (x)dx = J [q L (x) + q L (x)]dx. (13) 



The contributions of the Dirac continuum and of the discrete valence level to this moment are 
respectively 0.195 and 0.130, which shows that the former dominates over the latter. The sum 
of these two numbers gives L q ~ 0.325, which means that about 65 % of the nucleon spin is 
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Figure 1: (a) The theoretical predictions of the CQSM for the quark and antiquark orbital 
angular momentum distribution functions qi{x) and (b) the isosinglet quark polarization 
Au(x) + Ad(x). The long-dashed and dash-dotted curves respectively stand for the con- 
tributions of the discrete valence level and that of the Dirac continuum in the self-consistent 
hedgehog background, whereas their sums are shown by the solid curves. The distributions 
with negative x are to be interpreted as the antiquark distributions. 

carried by the orbital angular momentum of quarks and antiquarks at the energy scale of the 
model 



We can also know the separate contributions of quarks and antiquarks to the total orbital 
angular momentum. They are shown in table 1 together with the corresponding separation 
for the quark spin fraction of the nucleon, or the flavor-singlet quark polarization AS. One 
sees that the flavor-singlet polarization of antiquark is negative but fairly small in magnitude. 
(Note however that this does not necessarily mean that the polarizations of all the antiquarks 
are small. It has been shown that the CQSM predicts fairly large isospin asymmetry of 
the longitudinally polarized antiquark distribution in the nucleon, i.e. Ad{x) — Au(x) < 0.) 
On the contrary, a sizable amount (~ 30 %) of the total orbital angular momentum is seen to 
be carried by antiquarks with soft (small x) component. 

At this point, we want to make a short comment on the recent lattice QCD calculation of 



quark orbital angular momentum Aside from the fact that it is based on the so-called 
quenched approximation, one should also notice the fact that this is an indirect calculation 
of the quark orbital angular momentum L q as well as the total gluon angular momentum J g . 
They first extract the total quark angular momentum J q from the calculation of the quark 
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Table 1: The separate contributions of quarks and antiquarks to the first moment AE and L q 
at the scale of the model. 





quark 


antiquark 


total 


AS 


0.397 


- 0.047 


0.350 




0.229 


0.096 


0.325 


\ AS + L q 


0.427 


0.073 


0.500 



energy-momentum tensor form factors. Further combining the previous lattice calculation |30 
of the quark spin content AS, they obtain L q , thereby extracting J g from the total nucleon 
spin sum rule. This dictates that the extracted values of L q and J g would be sensitive to the 
uncertainties of J q and AS obtained in the numerical simulation. The main conclusion of this 
analysis is that about 25% of the nucleon spin originates from the quark intrinsic spin, while 
about 35% comes from the quark orbital angular momentum, which in turn means that the 
remaining 40% of the nucleon spin is due to the glue. If this is confirmed by more precise 
calculation in the future, the present investigation done with neglect of the explicit gluonic 
degrees of freedom needs a considerable amendment. At least, some of the quark contributions 
obtained within the QCSM must be redistributed to either or both of the gluon spin and orbital 
angular momentum. It is interesting to see that scaling the CQSM predictions (5 AE) ~ 35% 
and (L q ) ~ 65% by the factor of 0.6, we would obtain (| AS) ~ 21% and (L q ) ~ 39%, which 
is rather close to the corresponding numbers 25% and 35% obtained in the above lattice QCD 
study. An interesting common feature in both theories is the dominance of the quark orbital 
angular momentum over the intrinsic quark spin fraction. 

Now we are in a position to investigate the scale dependence of the full nucleon spin contents 
together with the corresponding distribution functions. As already mentioned, we assume that 
the gluon polarization and the gluon orbital angular momentum are both zero at the starting 
energy scale of Q 2 = Q 2 ni , i.e. 

AE(QL) = 0.350, Ag(Q 2 ni ) = 0, L q (Q 2 ini ) = 0.325, L g (Q 2 m ) = 0, (14) 



with the normalization 



^AE(Q 2 ) + Ag(Q 2 ) + L q (Q 2 ) + L g (Q 2 ) = \ . (15) 



) 2 

- in 

structure functions gi(x, Q 2 ), g™(x, Q 2 ) and gi(x,Q 2 ) based on the evolution equation at the 
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The question here is how to choose Qj ni - In our previous analysis of the longitudinally polarized 

J>f„ r>2\ „nt„ na\ n „A „d(„ r)2\ 



next-to-leading-order (NLO), a good agreement with the available high-energy data has been 
obtained with the choice Qf ni = 0.30 GeV 2 [^9| ,[pD[. Unfortunately, the evolution equation 
including the orbital angular momentum distributions is known only at the leading order (LO) 
Expecting that the use of a little smaller value of Q 2 ni would compensate the defect 
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of using this lower-order evolution equation, we proceed as follows. That is, we first solve 
the evolution equation for A£(Q 2 ) and Ag(Q 2 ) at the NLO ( in the MS scheme) by setting 
Qtni = 0.30 GeV 2 P^,PB|. Next, we solve the evolution equation for the full nucleon spin 
contents at the LO with use of a little smaller Q 2 inil i.e. Q 2 ni = 0.23 GeV 2 , which is determined 
so as to roughly reproduce the Q 2 -evolution of the gluon polarization Ag(Q 2 ) obtained by 
solving the NLO evolution with Q 2 ni = 0.30 GeV 2 . 



(a) 



0.8 



0.6 - 



0.4 



0.2 i 



Q 2 ini = 0.30 GeV 2 Ag(Q 2 ) 



AL(Q/) 



SMC data for A£(Q ) 



20 40 
Q 2 (GeV 2 ) 



60 




Q 2 (GeV 2 ) 



Figure 2: (a) The NLO evolutions of the flavor-singlet quark and gluon polarizations. The 
experimental data given at Q 2 = 10 GeV 2 corresponds to the NLO analysis by the SMC group 
M. (b) The LO evolution of the full nucleon spin contents. 



We first show in Fig.2(a) the NLO evolution for the quark and gluon polarization in com- 
parison with the latest SMC data for A£(Q 2 ). One can say that the prediction of the CQSM 
for AE is qualitatively consistent with the recent NLO analysis of the SMC group given at 
Q 2 = 10 GeV 2 §. One also sees that the gluon polarization rapidly grows with increasing Q , 
even if we have assumed Ag = at the initial energy scale of Q 2 ni = 0.30 GeV 2 . 

Shown in Fig.2(b) is the LO evolution of the full nucleon spin contents including the orbital 



angular momentum. They are calculated by using the known analytical solutions [31 



AS(Q 2 ) = AS(Q 2 ni ) = constant, 
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(16) 



MQ 2 ) = -33^^) + l (A,(QL) + 3 3^AE(g L ) 1, (17) 
W) = -^AS(QL) + ' 3 '"' 



2 ™" 2 16 + 3 ri/ 

2(16+3ny) 

+ ^-^ji / + 1 , 1 S„, 



to/ \ 2 2 16 + 3n/ 



1 1 6 

L g (Q 2 ) = -AgiCfiJ + - 



2 16 + 372/ 

2(16+3ny) 

£ \ "3(33-2^) ^ r „ ^ i a ^ 2 ^ 1 16 



t j v i.(«-) + AsKK-) - 5 ■ (19) 

with t = hiQ 2 / Aq CD , t = lnQ^/Ag^ and n/ = 3. The starting energy of this evolution 
is taken to be Q 2 ni = 0.23 GeV 2 , as mentioned above. As is widely-known, the flavor- singlet 
quark polarization AE is scale independent at the LO. Also relatively established is that the 
gluon polarization Ag(Q 2 ) shows a logarithmic growth with increasing Q 2 , while it is largely 
compensated by the similar decrease of the gluon orbital angular momentum L q (Q 2 ). In fact, 
it was pointed out by Ji et al. |3l|] that, as Q 2 — > oo, the sum of these two moments approaches 
its asymptotic value given by 

1 16 

J g = Ag + L a — ► , (20) 

q y 9 2 16 + 3n/ v ; 

which indicates (with ri/ = 3) that about 64% of the nucleon spin is carried by the gluon fields 

in the asymptotic energy. 

On the other hand, the evolution of quark orbital angular momentum is strongly dependent 

on the initial condition given at the low renormalization point. Since the scale-dependent factor 

2 (16+3 ni) 

^-3(33-2^ (2i) 

is a decreasing function of Q 2 (for 33 — 2n/ > 0), we conclude that L q (Q 2 ) is a decreasing 
function of Q 2 as far as 

L q (QL) + ^AE(QL) " 9 T^— ' ( 22 ) 
2 2 16 + 3 71/ 

is positive, while it is a increasing function of Q 2 if the same quantity above is negative. In 
the case of the CQSM, it gives L q (Q 2 ni ) + \ AE(Q 2 ni ) - \ = 0.325 + 0.175 - 0.18 > 

at Q 2 = Qj ni , so that it is consistent with the fact that L q (Q 2 ) shown in Fig.3 is a decreasing 
function of Q 2 . The Q 2 dependence of L q (Q 2 ) looks strong only in the relatively lower Q 2 
region, and beyond Q 2 ~ 5 GeV 2 it slowly approaches its asymptotic value given by 

L t[ Q>) _ -lAE(QL) + (23) 
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Substituting our model predictions AE(Qf ni ) = 0.35 with rif = 3, this gives L q — > 0.005, 
which denotes that the quarks barely carry orbital angular momentum at the very high energy 
scale. Note, however, that the precise fraction of the quark orbital angular momentum in the 
asymptotic domain depends on the delicate cancellation of the two terms in (23). 

It appears that the observation above contradicts one of the conclusions drawn in the recent 
paper by Scopetta and Vento [Q . They concluded that the quark orbital angular momentum 
can be important at large Q 2 , due to evolution, if it is not negligible at the scale of the hadronic 
model. It seems that this conclusion is mainly drawn from Fig. 5(b) in ||34|| , which shows that 
L q (Q 2 ) is an increasing function of Q 2 . This increasing behavior of L q (Q 2 ) appears to be 
inconsistent with the Q 2 -dependence of the quark orbital angular momentum distributions 
illustrated in Fig.4(a) of the same paper. In fact, the initial conditions of their "D-model" 
are given by AE(Q 2 ni ) ~ 0.4, Ag{Q 2 ni ) ~ 0.1, L q (Q 2 m ) ~ 0.145, and L g (Q 2 ni ) ~ 0.055. 
Since this gives L q (Q 2 im ) + ~ A£(Qf ni ) — \ ^+( n — 0.165 > 0, L q (Q 2 ) must be a decreasing 
function of Q 2 if our argument above is correct. Probably, the cause of discrepancy resides 
m an error m eqs. (17), (18) of @. The Q 2 -dependent factor &- 38 / 81 (with the definition 
b = a s (Q 2 )/as(Qj ni )) in (17) and (18) should be replaced by 6+ 50 / 81 for n f = 3. (Note that 
^-38/8i j g an i ncreas i n g function of Q 2 , while 6+ 50 / 81 is a decreasing function.) 

Next, we show in Fig. 3 the theoretical prediction of the CQSM for the spin and orbital 
angular momentum distribution functions and their LO evolutions. Here, we have used the 
LO evolution code provided by Martin et al. [p7fl . Here, the solid and long-dashed curves 
respectively represent x AS(i, Q 2 ) and x L q (x, Q 2 ), whereas the dashed and dash-dotted curves 
stand for x Ag(x,Q 2 ) and x L g (x,Q 2 ). At the initial energy scale of Q 2 ni = 0.23 GeV, we set 
Ag(x) = L g (x) = 0. On the other hand, AE(x) and L q (x) at the same scale are obtained as 

AE(i) = [Au(x) + Ad(x)] + [Au(x) + Ad(x)] , (24) 
L q (x) = q L (x) + q L (x), (25) 

from the predictions of the CQSM shown in Fig.l. The logarithmically increasing behavior 
of Ag(x, Q 2 ) and the logarithmically decreasing behavior of L g (x, Q 2 ) are just what has been 



pointed out in several previous papers p7|,p4"|. 



As pointed out in [p7|, L g (x,Q ) at high energy scale is rather insensitive to the input distri- 
butions for orbital angular momentum at the starting energy, and it is basically determined 
by the polarized quark singlet and gluon distributions. On the other hand, the behavior of 
L q (x, Q 2 ) at the moderate value of Q 2 is strongly dependent on the input orbital angular mo- 
mentum distribution. In fact, Fig.4 shows that L q (x,Q 2 ) remains positive at all the energy 
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Figure 3: The LO evolutions of the spin and orbital angular momentum distribution functions 
for quarks and gluons. the solid and long-dashed curves respectively represent x AE(x, Q 2 ) and 
x L q (x, Q 2 ), whereas the dashed and dash-dotted curves stand for x Ag(x, Q 2 ) and x L g (x, Q 2 ). 
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scales, although its magnitude gradually decreases as Q 2 increases. We point out that this be- 
havior of L q (x,Q 2 ) is consistent with the scale dependence of the corresponding first moment 
L q (Q 2 ) shown in Fig. 2(b). 

Clearly, the origin of this unique behavior of L q (x,Q 2 ) and L q (Q 2 ) can be traced back to 
the extraordinary large and positive orbital angular momentum for quarks and antiquarks at 
the hadronic scale predicted by the CQSM. The results above may be compared with several 
previous investigations of the scale dependence of the orbital angular momentum distributions, 
the input distributions of which are prepared in a somewhat arbitrary way. For instance, two 
different scenario have been investigated in p7[ . The first scenario called the GRSV standard 
scenario has the property that ~ AS(Qi„j) + Ag(Q 2 ni ) ~ 0.475 with very small orbital angular 
momentum for quarks and gluons, i.e. L q (Q 2 ni ) +L g (Q 2 ni ) ~ 0.015. The second rather extreme 
scenario assumes very large gluon polarization such that Ag(x, Qi n i) = g(x, Q 2 ni ), which in turn 
dictates large and negative orbital angular momentum such that L q (Q 2 ni ) +L g (Q 2 ni ) = — 0.43. 
A common feature of these scenarios is that the evolved L q (x,Q 2 ni ) has large and negative 
values in wide range of x, which is quite different from the above-mentioned result of the CQSM 
shown in Fig. 4. Unfortunately, any of these quite different scenarios cannot be excluded on 
the basis of our present knowledge, since all of these reproduces the available data for the 
longitudinally polarized deep-inelastic structure functions for the nucleon and the deuteron 
at least qualitatively. Highly desirable is some direct experimental information for any of 
Ag(x,Q 2 ), L q (x,Q 2 ) and L g (x,Q 2 ). 



4 Conclusions 

An outstanding feature of the CQSM is that it can explain all the qualitatively noticeable 
features of the recent high-energy deep-inelastic-scattering observables including NMC data 
for F^{x) — F£(x), F^{x) / F^x) |36| , the Hermes |35| and NuSea |37j data for d(x) —u(x), the 



EMC and SMC data for gi(x), and gf(x) [2T|,0], without any free parameter except for 



the starting energy scale of the DGLAP evolution equation p9[,P0||. In the present paper, we 
have evaluated the spin and orbital angular momentum distribution functions within the same 
theoretical framework, and have investigated their scale dependences. It has been shown that 
the model can explain the smallness of the quark spin fraction at least qualitatively without 
assuming large gluon polarization at the low renormalization point. A key ingredient here is 
the large and positive orbital angular momentum carried by quarks and antiquarks with soft 
(small x) component. The predicted large and positive quark orbital angular momentum at 
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the low renormalization point necessarily affects the scenario of the scale dependence of the 
nucleon spin contents. The solution of the LO evolution equation indicates that the quark 
orbital angular momentum is a decreasing function of Q 2 and it may play insignificant role at 
the asymptotic energy scale, although the rate of reduction is quite slow beyond Q 2 ~ 5 GeV 2 . 
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